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摘要 ”特征 线 方法 (Method of Characteristics; MOC) 因 其 具备 强大 的 几何 处 理 能 力 ， 且 在 计算 过 程 中 亦 能 兼顾 
计算 成 本 和 计算 精度 ， 被 广泛 应 用 于 高 保 真 数值 模拟 计算 中 。 常 见 的 输 运 计算 方法 除 MOC 外 ,还 包括 碰撞 概率 
ik (Collision Probability Method, CP) 和 界面 流 法 〈Interface Current Method, IC) 等 。 本 文 从 方法 理论 以 及 数 
值 计算 上 将 MOC, CP 和 IC 进行 比较 分 析 ， 对 于 评估 其 在 pin-by-pin 计算 中 的 能 力 尤 为 重要 。 同 时 在 MOC it 
算 中 , 不 同 的 参数 选择 会 对 计算 成 本 和 计算 精度 产生 影响 , 因此 有 必要 进行 敏感 性 分 析 以 寻求 最 佳 参数 。 本 文江 
先 将 三 种 计算 方法 从 原理 上 进行 比较 分 析 ， 再 基于 2D C5G7-MOX 基准 题 完 成 了 数值 计算 及 MOC 参数 敏感 性 
初步 分 析 。 计 算 结 果 表 明 : MOC 在 计算 精度 、 计 算 效 率 和 内 存 开 销 上 均 优 于 CP 和 IC MOC 的 计算 耗 时 和 内 
存 开销 分 别 为 23.9min 和 37.5MB， 与 参考 解 的 相对 误差 仅 为 6.04x104。 而 CP 和 IC 的 计算 耗 时 分 别 为 MOC 的 
56.7 倍 和 15.6 倍 ， 内 存 开销 分 别 为 MOC 的 407.7 倍 和 32.8 倍 。 进 一 步 通过 参数 敏感 性 分 析 发 现 ， 网 格 划 分 对 
计算 内 存 开 销 以 及 计算 时 间 的 影响 最 大 , 极 角 的 选择 对 计算 精度 的 影响 最 大 , 并且 给 出 一 组 建议 参数 : 网 格 划分 
6x6， 极 角 为 GAUS 且 数 目 为 2， 方位 角 个 数 为 30。 该 组 参数 的 计算 耗 时 为 4$5.4min， 内 存 开销 为 264.7MB， 相 
对 误差 为 $.9x105， 归 一 化 后 的 栅 元 均 方 根 误 差 为 0.00255 。 
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Abstract [Background]: Method of Characteristics (MOC) is widely used in high-fidelity numerical simulation 
because of its powerful geometric processing ability and its ability to take into account calculation cost and accuracy. 


In addition to MOC, the common transport calculation methods also include Collision Probability Method (CP) and 
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Interface Current Method (IC). [Purpose]: In this paper, the comparative analysis of MOC, CP and IC from method 
theory and numerical calculation is particularly important for evaluating their ability in pin-by-pin calculation. At the 
same time, in MOC calculation, different parameter selection will affect the calculation cost and calculation accuracy, 
so it is necessary to carry out sensitivity analysis to find the best parameter. [Methods]: In this paper, the three 
calculation methods are compared from the principle, and then the numerical calculation and the preliminary analysis 
of the sensitivity of MOC parameters are completed based on the 2D C5G7-MOX reference problem. [Results]: The 
calculation time and memory cost of MOC are 23.9min and 37.5MB respectively, and the relative error between MOC 
and reference solution is only 6.04x10^*. The computing time of CP and IC is 56.7 times and 15.6 times that of MOC, 
and the memory cost is 407.7 times and 32.8 times that of MOC, respectively. Based on the sensitivity analysis of MOC 
parameters, a set of suggested parameters are given: grid division is 6x6, pole Angle is GAUS and number is 2, and 
azimuth Angle is 30. The calculation time of this set of parameters is 45.4min, the memory cost is 264.7MB, the relative 
error is 5.9x10?, and the normalized RMS error is 0.00255. [Conclusions]: The results show that MOC is superior to 
CP and IC in accuracy, efficiency and memory cost. Through the parameter sensitivity analysis, it is found that the grid 
division has the greatest influence on the calculation memory cost and calculation time, and the choice of polar Angle 


has the greatest influence on the calculation accuracy. 
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随 着 新 型 反应 堆 的 发 展 ， 更 多 非 均匀 性 强 
来 了 挑战 。 在 高 保 真 反应 堆 物 理 计 算 中 ， 常 ) 


储 型 出现， 为 堆 芯 物理 计算 与 分 析 带 
包括 : 特征 线 方法 (Method of 


Characteristics: MOC) ~, fil fit EAS 7 C Collision Probability Method, CP) 和 界面 流 法 (Interface Current Method, 


IC) 等。CP 和 IC 方法 都 属于 积分 输 运 方法 ，i 


算 栅 元 的 热 : 


反应 堆 物理 计算 中 的 一 种 引 
Bl。 当 下 积分 输 运 方法 作为 中 子 输 i 


5B5l, MOSRA!6AE , 


子 利 用 系数 023， 并 且 得 到 了 很 高 的 计算 精度 。 因 


日 A.Amouyal、P.Benoist 和 J.Horowitz M.F] T il 
此 积分 输 运 方法 曾 在 相当 一 段 时 间 内 被 作为 
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HfFBeuE E 
以 实现 ， 其 中 包括 APOLLO, WIMSD- 


要 方法 ， 特 别 是 如 


特征 线 方法 最 早 由 英 


升 ， 特 征 线 方法 凭借 其 精确 的 数值 计生 
1。 目 前 国内 外 许多 程序 都 采 


法 044， 该 方法 巧妙 地 运 月 


国学 者 AskewD 应 用 于 求解 复杂 的 输 


被 广泛 应 用 


运 计算 ， 随 着 近 些 年 来 计算 机 性 能 的 不 断 攀 


和 灵活 的 儿 何 特性 ， 在 高 保 真 反应 堆 物 理 分 析 中 得 到 了 广泛 的 应 用 


解 03。 特 征 线 方法 具备 强大 的 几何 处 理 能 力 ， 能 将 角度 变量 和 空间 变量 
几何 ， 并 且 由 于 特征 线 中 线 与 线 之 间 的 独立 性 ， 使 得 特征 线 方 ; 
顾 计算 效率 和 计算 精度 。 MOC、CP、IC 三 种 方法 各 具 
的 挑战 ， 因 此 有 必要 对 以 上 三 种 输 运 计算 方法 进行 分 析 比 较 ， 评 佑 其 


征 线 方法 作为 输 运 求解 器 ， 包 括 DRAGON[910 
等 。 特 征 线 方法 (Method of Characteristicss MOC) 是 沿 中 子 轨 迹 以 特征 线形 式 求解 中 子 输 运 方程 的 一 种 方 


、CASMOULI2、WIMSI03] 


昌 了 数学 技巧 ， 将 中 子 输 运 方程 


微分 形式 表示 ， 沿 着 中 子 运行 的 轨迹 进行 积分 求 


司 时 进行 处 理 ， 理 论 上 能 处 理 任 意 
备 天 然 的 并 行 能 力 ， 在 计算 过 程 中 亦 能 


真 计 算 时 ， 都 会 面临 计算 成 本 
在 pin-by-pin 计算 中 的 能 


本 文 首先 比较 了 特征 线 方法 (Method of Characteristics: MOC )、 碰 撞 概 率 法 (Collision Probability Method, 


CP) 和 界面 流 法 〈Interface Current Method, IC) 的 计算 原理 ， 并 基于 2D C5G7-MOX 基准 题 分 析 了 三 种 方 


法 的 计算 精度 、 


WU. Wf 
极 角 、 方 位 角 以 及 网 格 划分 3 


A 


耗 时 、 计 算 精 度 和 内 存 3 


计算 效率 和 内 存 3 


题 的 参数 选择 


与 优化 提供 


征 线 方法 本 身 在 计算 的 过 程 中 ， 不 同 参数 的 组 合 会 导致 
效率 和 内 存 开销 发 生变 化 ， 需 要 进一步 开展 MOC 参数 敏感 性 分 析 。 
于 展 了 参数 敏感 性 分 析 ， 但 仅 丰 
生 分 析 中 的 参数 类 型 ， 对 部 分 参数 的 选取 范 
开销 ， 能 够 较为 全 


中 Miral 6 A EJ 


究 了 上 述 三 个 参数 对 计算 精度 的 影响 。 本 文 完 
进行 了 深入 探究 ， 分 析 比 较 了 各 种 参数 的 计算 
的 为 特征 线 方法 在 计算 2D C5G7-MOX 这 类 强 非 均匀 性 组 件 问 


ERME: 基于 C5G7-MOX 的 中 子 输 运 计算 方法 比较 及 MOC 参数 敏感 性 初步 


分 析 


1 计算 方法 


1.1 特征 线 方法 


特征 线 方法 (Method of Characteristics, MOC) 的 核心 思想 是 沿 着 特征 线 方向 进行 积分 ， 将 中 子 输 运 方 


程 转化 成 特征 线 方程 ， 再 通过 对 每 一 条 平行 特征 线 进行 输 运 求解 ， 从 而 完成 对 整个 目标 区 域 的 计算 。 从 理 


论 上 讲 ， 特 征 线 方法 具有 优异 的 几何 适应 能 力 ， 只 要 特征 线 足够 密 ， 便 可 对 任意 形状 的 几何 区 域 进行 输 运 


求解 。 


使 用 特征 线 方法 进行 中 子 输 运 方程 求解 时 ， 需 将 待 求解 问题 的 几何 区 域 划分 成 众多 平 源 区 ， 并 假设 在 
所 划分 的 平 源 区 中 ， 宏 观 截 面 和 中 子 源 项 均 为 常数 。 如 图 1 所 示 ， 对 于 任意 平 源 区 六， 其 中 带 箭 头 的 平行 实 
线 分 别 为 穿 过 该 区 域 的 特征 线 ， 两 虚线 之 间 的 距离 为 第 条 特征 线 的 宽度 记 作 564mnxk， 特 征 线 沿 厦 方位 角 


为 m、 极 角 为 nx"， 其 中 截面 和 源 项 近似 分 别 如 式 CIO 和 式 〈2) 所 示 。 
YXus) = Xu 
Qim(s) 入 Qim 


(1) 
(2) 


nea uic du a Led 
化 的 源 项 ， 于 是 稳 态 多 群 中 子 输 运 方程 可 以 表示 为 品 


sinOn, 2 Qunnk (simnx) t NtiDisunk (simnx) Qim 


区 内 不 随 位 置 而 发 生变 


(3) 


AX O0 中 : 为 中 子 沿 着 Qmg 方 向 飞行 的 极 角 ;simmx 为 特征 线 在 x-y 平面 投影 所 得 的 长 度 ，@imnk(5) 为 


中 子 角 通 量 。 式 (3) 为 一 阶 线性 微分 方程 ， 可 以 解析 求解 得 到 下 式 : 


uk = Pin nexp(—Deisimnr) tan Lim mji mi exp(- Yeisimnk)] (4) 


X (4) 中 : p83tiix 为 出 口 处 的 中 子 角 通 d mE. XX (2) 沿 着 特征 


线段 进行 积分 ， 由 中 子 守恒 关系 可 得 : 


Pimnx Simn,k 7 pen $i n keXb(7XtiSimnk) TY Pun mji 7 exp(- Fiisimnx)|ads (5) 


人 


in out 
Qim m Pimnk— Pimn,k 


Pimnk = Yu 


XtiSimnk 


nw 


pny 
图 1 特征 线 方法 示意 图 


Fig.l Diagram of method of characteristics 


对 在 区 域 i 内 的 以 方位 角 m、 极 角 n 的 所 有 特征 线段 ， 以 体积 为 权重 进行 求 和 ， 
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(6) 


得 其 到 区 域 i 内 的 平均 
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角 通 量 密度 ; 
Bimn = Penna 5 
R CD 中 : dua EN m ifc. m BUG SERE i KRA RERBA, VRR DOR ER 


大 小 ， 对 其 方位 角 和 极 角 进行 加 权 平 均 ， 可 得 i 区 域内 的 中 子 通 量 为 : 
pi = 5a Om DN On Pimn (8) 
式 (8) P: wm 表示 方位 角 求 积 系数 ，M 为 方位 角 的 离散 数目 ，w 为 极 角 求 积 系 数 , N 为 极 角 的 离散 数目 ， 
将 式 CD 和 式 (6) RAR (8) 中 可 得 ，i 区 域内 的 标量 中 子 通 量 密度 如 下 : 
angi | EM- Om Xni On Óima Ota PAmnk 
$i T A + Ec es nl T.I (9) 
根据 式 (9) 可 知 ， 在 计算 各 平 源 区 的 中 子 通 量 密度 时 ， 极 角 n、 方 位 角 m、 特 征 线 密度 特征 线 宽度 
54mmxk 以 及 各 平 源 区 域 划分 的 精细 程度 都 会 直接 影响 中 子 通 量 的 计算 精度 。 因 此 ， 本 文 基 于 上 述 参数 开 
展 敏感 性 分 析 ， 比 较 不 同 参数 组 合 对 计算 精度 、 计 算 时 间 和 内 存 开 销 的 影响 程度 ， 为 后 续 特 征 线 参 数 的 选 
择 和 优化 提供 参考 。 


1.2 碰撞 概率 法 

碰撞 概率 法 〈Collision Probability Method, CP) 是 求解 积分 中 子 输 运 方程 的 最 基本 方法 ， 运 用 碰撞 概 
率 法 进行 求解 时 ， 将 系统 划分 成 个 互 不 相交 的 均匀 子 区 ，V=1,2,……N， 其 中 当 ij ViVeo. BT 
Vi 划分 的 足够 小 ， 可 以 近似 地 认为 在 每 一 子 区 域内 的 截面 参数 等 于 常数 。Pi; 为 j 区 内 产生 的 一 个 各 向 同性 
中 子 ， 不 经 任何 碰撞 到 达 i 区 发 生 首 次 碰撞 的 概率 L717， 由 于 Pj 的 求解 与 j 区 的 源 中 子 或 中 子 通 量 密度 的 分 
布 有 关 ， 而 中 子 通 量 密度 正 是 所 要 求 的 未 知 量 。 因 此 需要 引入 平 源 近 似 ， 即 当 子 区 域 划分 得 较 小 时 ， 近 似 
也 认为 每 个 区 域 的 源 中 子 或 中 子 通 量 密度 的 分 布 等 于 常数 ， 同 时 为 了 简化 计算 ， 亦 要 假设 入 射 中 子 在 子 区 
表面 上 的 空间 分 布 是 均匀 的 ， 且 入 射 中 子 的 角 分 布 各 向 同性 09。 故 在 多 群 近似 的 条 件 下 ， 便 可 对 局 ) 以 及 
通 量 进行 求解 如 下 : 


x 


PLE efre -7)] i 
Pus mp de A T a Y (10) 
— 51I G E 
BeigbigV = Xj [g's1 ggpig' + Sio|PuaV; ,9 = 12G (11) 


X (10) 中 : Pijg 为 g 群 内 在 j 区 内 产生 的 各 向 同性 的 中 子 ， 不 经 任何 碰撞 到 达 i 区 发 生 首次 碰撞 的 概率 。 
Fig 为 g 群 中 子 鹤 面 近似 得 到 的 宏观 截面 ，r” 一 + 为 空间 上 位 置 的 变化 ，tg (r” 一 + ) 为 从 r' 到 7 点 之 间 的 
自由 程 ，V 为 区 域 j 的 体积 。 式 (11) 称 为 碰撞 概率 形式 的 中 子 积分 输 运 方程 ， 其 中 等 号 右 端 方 括号 内 项 表 
示 在 态 内 单位 体积 内 g 群 的 源 中 子 数 ， 因 此 ， 等 号 右 端 整体 便 表示 所 有 7 区 的 g 群 中 子 在 Vi 区 域内 首次 发 生 
碰撞 的 数目 ， 等 于 等 号 左 端 丰 内 的 碰撞 数 品 ;ygis 灰 。 且 在 整个 碰撞 概率 法 求解 的 过 程 中 ， 中 子 飞 行 首次 磁 
撞 概 率 和 矩阵 忆 的 求解 最 为 复杂 和 耗 时 品 。 


1.3 界面 流 法 

界面 流 法 (Interface Current Method, IC) 的 基本 思想 是 : 将 系统 分 为 若干 子 区 ， 对 每 个 子 区 进行 均匀 
化 处 理 后 ， 再 对 每 个 子 区 建立 中 子 输 运 方程 ， 各 子 区 之 间 通 过 界面 上 的 中 子 流连 续 的 条 件 进 行 耦合 ， 以 此 
来 求解 中 子 通 量 。 

将 系统 划分 为 NIXNJ 个 子 区 ， 假 定 中 子 源 的 角 分 布 为 各 向 同性 ， 空 间 分 布 采 用 线性 分 布 ， 即 所 谓 的 线 
性 源 近似 ， 于 是 中 子 源 项 可 表示 为 07; 


y M 2 G| +a (2-2) +a (2-2, 0x «0v «o (12) 


XX C2) P: Q(x,y,0) 为 子 区 内 点 (x,y) 人 处 沿 0 方 向 的 中 子 源 密度 ，Q 为 子 区 的 平 源 强 ，a、b 为 子 区 在 x、y 
JARAK, ap ay HME x y 方向 上 的 线性 分 布 系数 。 
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如 图 2 的 矩形 几何 为 例 ， 对 于 子 区 (i, 让 设 体积 为 Vi;， 根 据 中 子平 衡 可 得 在 界面 流 方程 为 : 
ZKk=4BCD zz Oa Jia — X42 Jia) + «goij = Qu (13) 


X 030 P: 9ij 为 子 区 G, D RFID, LATE G D 的 总 截面 ，Qi T IX G, D 内 的 平均 源 
强 ( 包 括 散 射 源 、 和 裂变 源 、 外 中 子 源 )。E 为 矩形 几何 子 区 的 四 个 界面 CA, B, C, DO 如 图 2 所 示 ，9 为 各 
界面 上 中 子 飞 行 方向 划分 为 四 个 象限 ， 每 个 边界 的 四 个 象限 在 图 2 中 均 己 标 出 ， 厌 ,为 子 区 界面 上 象限 q 的 
出 射 中 子 流 ， Jea FX AH k 象限 q 的 入 射 中 子 流 。 

界面 流 法 在 求解 通 量 时 ， 考 虑 的 是 每 个 子 区 域内 源 与 相 邻 子 区 域 的 边界 的 中 子 守恒 关系 271， 每 个 子 区 
只 需 计 算 首 次 穿 透 和 泄露 两 个 概率 ， 且 计算 的 概率 数目 只 与 系统 内 所 出 现 的 子 区 类 型 成 正比 。 而 碰撞 概率 
法 中 首次 磁 撞 的 计算 量 ， 却 是 与 系统 内 子 区 数目 的 平方 成 正比 。 因 此 界面 流 法 相 较 于 碰撞 概率 法 ， 具 有 较 
高 的 计算 效率 。 


LI 


图 2 象限 的 划分 与 界面 标号 


Fig.2 Quadrant division and interface labeling 


2 计算 模型 


2.1 2D C5G7 MOX 基准 题 
Reflection Boundary 


Reflection Boundary 
Vacuum Boundary 


Vacuum Boundary 
3 C5G7-MOX 组 件 布置 示意 图 
Fig.3 C5G7-MOX assembly layout diagram 


C5G7 基准 题 是 由 OECD/NEA 于 2003 年 发 布 的 基准 题 ， 主 要 用 于 评价 现 有 输 运 计算 方法 进行 全 堆 芯 
非 均 匀 计 算 的 能 力 09。 该 基准 题 包 含 y UOM MOX 两 种 燃料 组 件 ， 每 个 组 件 的 边 长 为 21.42cm， 由 17x17 
的 燃料 栅 元 组 成 。 外 围 慢 化 剂 区 域 的 宽度 与 组 件 边 长 相同 。2D C5G7-MOX 基准 题 的 结构 如 图 3 Hz. dh 
侧 和 下 方 为 真空 边界 ， 左 侧 和 上 方 为 全 反射 边界 19。 各 类 组 件 中 的 各 个 栅 元 都 具有 相同 的 几何 参数 。 栅 元 
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几何 结构 如 图 4 所 示 , 其 中 包括 燃料 和 包 壳 打 混 的 圆 形 区 域 以 及 外 部 慢 化 剂 的 矩形 区 域 , 栅 元 边 长 为 1.26cm， 
燃料 区 域 的 半径 为 0.54cm。 

C5G7-MOX 基准 题 共 由 七 种 材料 构成 : UOo 燃料、4.3%、7.0% 和 8.7% 的 MOX 燃料 、 裂 变 室 、 导 向 管 
和 水 。 其 中 水 不 仅仅 用 于 填充 外 部 慢 化 剂 区 域 ， 还 用 于 填充 每 个 机 元 内 部 棒 外 的 区 域 ,， UOs 和 MOX 组 件 
内 的 燃料 棒 布置 如 图 5 所 示 。 


Moderator 


Fuel-Clad Mix 


1.08 
1.26 


4 C5G7-MOX 栅 元 结构 示意 图 
Fig.4 C5G7-MOX cell structure diagram 


图 5 C5G7-MOX 燃料 棒 布 置 示意 图 
Fig.S C5G7-MOX fuel rod layout diagram 


22 DRAGONS.0.8 

DRAGON 是 由 加 拿 大 蒙特 利 尔 理工 大 学 核 工程 研究 所 开发 的 , 基于 确定 论 方法 的 反应 堆 栅 格 计算 程序 
外 。 该 程序 采用 特征 线 方法 、 界 面 流 法 、 碰 撞 概 率 法 等 多 种 数值 方法 求解 中 子 输 运 方程 PZ0， 已 被 广泛 用 于 反 
应 堆 物 理 计 算 研 究 中 。DRAGON 具有 较 强 的 几何 处 理 能 力 ， 能 处 理 常 见 的 二 维和 三 维 几何 模型 。 基 于 
DRAGON 计算 2D C5G7-MOX 基准 题 的 方法 如 图 6 所 示 。 首 先 通过 MAC 模块 读 入 该 基准 题 中 各 材料 的 7 
群 截面 参数 ， 几 何 上 采用 51x51 的 栅 元 阵列 排 布 ， 其 中 顶部 和 左 侧 的 边界 条 件 为 镜面 边界 ， 右 侧 和 底部 的 
边界 条 件 为 真空 边界 〈GEO 模块 ) ， 然 后 对 待 求解 的 几何 区 域 进行 特征 线 布 置 与 追踪 ， 获 得 各 平 源 区 内 特 
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内 的 中 子 通 量 密度 及 反应 率 ， 完 成 二 维 输 运 计算 (ASM 模块 、FLU 模块 和 EDI 模块 ) 。 


Macrolib 
Geometry 


Ameliorated 
track 
Asmpij 


图 6 基于 DRAGON 的 特征 线 方法 
Fig.6 Method of characteristics based on DRAGON 


Geometry 
track 


Flux 


3 计算 结果 及 分 析 


本 文采 用 蒙特 卡 罗 程 序 OpenMC 的 计算 结果 作为 参考 解 。OpenMC 是 2011 年 由 美国 肪 省 理工 学 院 开 


征 线 的 相关 几何 信息 CNXT 模块 和 MCCGT 模块 ) ， 最 后 逐一 对 每 条 特征 线 进行 输 运 求解 ， 获 得 各 平 源 区 


TUN 


发 的 蒙特 卡 罗 ! 


子 和 光子 输 运 模拟 程序 ， 专 注 于 反应 堆 建 模 和 反应 堆 物 理 方法 研究 ， 旨 在 实现 高 保 真 


性 能 、 大 规模 仿真 计算 21223， 有 具备 临界 计算 、 固 定 源 计算 等 多 种 功能 23]。 
测试 环境 为 PC 单机 ， 有 具备 5.0GHz 十 四 核心 Intel Ef i9 处 理 器 ，16GB 内 存 。 通 过 OpenMC 共计 算 


^ [可 


模拟 800 代 ， 去 掉 前 200 代 ， 统 计 后 600 代 ， 每 代 跟 踪 600000 个 中 子 ， 以 得 到 各 燃料 栅 元 的 归 一 化 功率 及 


X (14) P: k 


D 为 本 文 计算 的 有 效 中 子 增殖 因子 ， 而 Ko 


有 效 中 子 增殖 因子 (kef 王 1.18652+0.00004)。 相 对 误差 计算 公式 如 下 : 


_ [kekse 
RE — e 4 e 


keff 


eff 


为 OpenMC 计算 的 有 效 中 子 增殖 因子 。 通 过 


(14) 
比较 


组 件 中 每 个 栅 元 的 功率 计算 误差 ， 从 而 进一步 对 功率 的 空间 分 布 进行 评估 ， 均 方 根 误差 计算 公式 如 下 : 
RMSE = [Eco 
N 


X (15) P: N 为 栅 元 的 总 数目 


的 归 一 化 功率 。 


3.1 不 同 中 子 输 运 方法 计算 比较 


首先 将 特 生 


E 线 方法 (MOC) 与 碰撞 概率 法 CCP) 、 界 面 流 法 AC) ， 通 过 对 2D C5G7-MOX 基准 


(15) 


，pp ;为 本 文 计算 所 得 第 ;个 栅 元 的 归 一 化 功率 ，po; 为 参考 解 中 第 ;个 栅 元 


表 1 展示 了 分 别 采用 特征 线 方法 、 碰 撞 概 率 法 和 界面 流 法 在 计算 2D C5G7-MOX 基准 题 的 耗 时 、 


开销 以 及 相对 误差 。 计 算 时 间 方 面 : 特征 线 方法 耗 时 最 短 ， 界 面 流 法 次 之 ， 碰 撞 概 率 法 最 长 。 其 中 界 
法 耗 时 是 特征 线 方法 的 15.6 倍 ， 碰 撞 概 率 法 耗 时 是 特征 线 方法 的 56.7 倍 ; 内 存 开销 方面 : 特征 线 方法 
存 开 销 远 小 于 碰撞 概率 法 和 界面 流 法 ， 其 中 碰撞 概率 法 内 存 开销 最 大 。 界 面 流 法 的 内 存 开销 是 特征 线 方法 
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题 计 
算 的 计算 精度 、 计 算 效 率 和 内 存 开 销 进 行 比较 。 再 从 三 种 方法 的 计算 原理 进行 解释 ， 能 够 更 全 面 地 对 特征 
线 方法 、 础 撞 概 率 法 、 界 面 流 法 的 计算 能 力 进 行 评 估 。 


内 存 
apo 
的 内 
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的 32.8 fi, RÉJERBEXSTERUVS 103 


E 


开销 是 特征 线 方法 的 407.7 倍 ; 


法 与 参考 解 相对 误差 高 达 3.8996x10-?。 


采用 碰撞 概率 法 计算 虽然 
若干 区 域 ， 每 个 区 域 的 中 子 通 量 密度 彼此 相互 耦合 ， 得 到 中 子 通 量 需要 求解 所 有 
这 些 都 会 使 计算 规模 变 大 ， 从 而 导致 耗 时 增多 内 存 3 


也 足够 精 丰 


但 是 会 消耗 大 量 时 


计算 精度 方面 特征 线 方法 的 计算 结果 最 为 
甫 确 ， 与 参考 解 的 相对 误差 只 有 6.04x104， 碰 撞 概 率 法 算得 的 结果 与 参考 解 相对 误差 6.14x104， 而 界面 流 


碰撞 概率 法 将 系统 分 成 


FRE PE, MAEI 


于 销 增 大 。 界 面 流 法 在 计 外 


区 域 的 首次 碰撞 概率 24。 


过 程 中 ， 由 于 中 子 流 过 各 子 


区 域 之 间 的 界面 必须 在 连续 条 件 下 耦合 ， 且 每 个 子 区 域 上 只 与 相 邻 子 区 域 耦合 ， 


只 需要 足够 密 的 特征 线 将 所 有 计 
此 可 以 通过 调整 特征 线 的 疏 密 程 度 ， 在 保证 精度 的 同时 也 能 减 4 
计算 2D C5G7-MOX 基准 题 这 类 强 非 均匀 性 的 问题 时 ， 计 算 精 度 、 训 


概率 法 和 界面 流 法 。 


-« 


x1 2D C5G7-MOX 燃料 组 件 计 算 结 果 


计算 方法 耗 时 

Computing method Elapsed time / min 

MOC 23.9 

CP 1355.6 

IC 374.1 
32MOC 参数 敏感 性 分 析 


进一步 对 影响 特征 线 方法 计生 


Tablel 2D C5G7-MOX fuel assembly calculation results 


内 存 开销 


Memory overhead / MB 


37.5 


15291.1 


1233.4 


因此 在 平 源 近 似 假设 下 计生 
强 非 均匀 性 问题 的 误差 较 大 2529, 会 使 其 误差 远大 于 特征 线 方法 和 碰撞 概率 法 ,特征 线 方 法 在 计算 过 程 中 ， 
算 区 域 进行 覆盖 ， 沿 着 特征 线 进行 积分 ， 从 而 完成 对 整个 区 域 的 计算 。 因 
计算 成 本 。 综 合 以 上 分 析 ， 特 征 线 方法 在 
| 算 效 率 和 内 存 开销 总 体 上 是 优 于 碰撞 


相对 误差 


Relative error / 10? 


60.4 
61.4 


3899.6 


结果 的 关键 参数 进行 敏感 性 分 析 。 由 式 〈9) 可 知 极 角 、 方 位 角 、 特 征 线 


密度 的 选择 都 会 对 中 子 通 量 的 计算 结果 产生 影响 ， 且 每 个 栅 元 划分 的 网 格 数目 也 会 对 计算 精度 、 计 算 效率 


和 内 存 开销 产生 影响 。 


R 2 Bin: 


参数 类 型 
Parameter type 
极 角 类 型 
Polar Angle type 
极 角 数目 
Number of polar angles 
方位 角 

Azimuth angles 

线 密度 


linear density / cm'! 


通过 改变 网 格 数目 、 方 位 角 、 极 角 、 特 征 线 密度 进行 计生 
数 对 计算 效率 以 及 精度 的 影响 从 大 到 小 为 : 网 格 数 
算 产 生 的 影响 大 小 顺序 ， 分 别 对 网 格 数 


， 共 得 到 44 组 数据 ， 初 步 分 析 得 出 ， 以 上 参 


目 、 极 角 、 方 位 角 、 特 征 线 密度 进行 计生 


表 2 参数 类 型 与 选择 区 间 


Table2 Parameter type and selection interval 


选择 区 间 


Selection interval 


GAUS, LCMD. CACB 


8. 12. 16. 30. 32. 64 


8. 10. 12. 15. 20 
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目 、 极 角 、 方 位 角 、 特 征 线 密度 。 因 此 ， 本 文 按照 对 计 


分 析 ， 参 数 类 型 及 区 间 如 


ERE: 基于 C5G7-MOX 的 中 子 输 运 计算 方法 比较 及 MOC 参数 敏感 性 初步 分 析 


网 格 划 分 


Mesh generation 


2x2. 4x4. 6x6. 8x8 


3.2.1 网 格 划 分 与 极 角 


400000 


200000 


Memory Overhead / kB 


Memory Overhead / kB 
N N N 
S 
S 


—a— CACB 2 "E, 
*— CACB 3 nad i 


4— GAUS 2 


4— LCMD3 


36 
v— GAUS 3 
*— LCMD2 

—4— LCMD 3 


16 36 64 
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图 7 不 同 网 格 和 极 角 下 的 内 存 开销 
Fig.7 Memory overhead for different grids and polar angles 


图 7 展示 


了 选择 不 同 的 网 格 划分 以 及 不 同 极 角 进行 计算 所 得 的 内 存 开销 ， 其 


The Number of Individual Cell Grids 


FERA es ET BIET 


格 数目 ， 纵 坐标 表示 整个 计算 过 程 得 内 存 开销 《单位 kB)， 其 中 不 同 的 曲线 分 别 代表 六 种 极 角 。 从 图 7 可 
知 ， 随 着 单个 栅 元 的 网 格 划 分 数目 的 增多 ， 内 存 开 销 也 随 之 显著 增 大 。 而 相 比 之 网 格 划分 对 内 存 开 销 的 影 
响 ， 不 同 极 角 的 选择 对 内 存 开销 的 影响 几乎 可 以 忽略 不 计 。 在 网 格 划分 为 2x2 的 


情况 下 ， 六 种 极 角 分 别 计 
算 的 内 存 开 销 为 3.8x10 和 kB; 在 网 格 划分 为 4x4 的 情况 下 ,六 种 极 角 分 别 计算 的 内 存 开 销 为 1.27x105kB; 在 


网 格 划 分 为 6x6 的 情况 下 ， 六 种 极 角 分 别 计算 的 内 存 开销 为 2.71x105kB;， 在 网 格 划 分 为 8x8 的 情况 下 ， 六 
种 极 角 分 别 计算 的 内 存 开销 为 4.36x105kB。 根 据 不 同 网 格 划 分 和 极 角 下 的 内 存 开 销 可 得 ， 相 同 网 格 数目 下 
的 不 同 极 角 划分 对 内 存 的 开销 影响 仅 在 400kB 以 内 ， 而 不 同 网 格 划分 对 内 存 开销 影响 非常 明显 ， 网 格 划 分 


越 密 内 存 开销 


更 越 大 ，8x8 网 格 划 分 比 2x2 网 格 划分 内 存 开销 多 3.98xl05kB， 是 
10000 
8000 上 
A 
s pA 
2 6000 上 Z 
E f 
$ 1000 F =— CACB 2 
a xA  |-4— CACB3 
- = —v— GAUS2 
20 r- P + GAUS 3 
t 4— LCMD2 
»— LCMD 3 
0 L fi 1 fi 
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图 8 不 同 网 格 和 极 角 下 的 总 耗 时 
Fig.8 Total time spent at different grids and polar angles 
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2x2 网 格 划 分 的 11.4 倍 。 


E 技 Ñ 


20XX, XX: XXXXXX 


图 8 展示 了 不 同 网 格 和 极 角 下 的 耗 时 ， 其 中 横 坐 标 表示 单个 栅 元 网 格 数目 ， 纵 坐标 表示 整个 计算 过 程 


的 耗 时 (单位 s), 其 中 不 同 的 线 分 别 代 表 六 种 极 角 。 从 图 


8 可 知 , 网 格 划 分 对 耗 时 的 影响 要 大 于 极 角 选择 ， 


且 随 着 网 格 划分 数目 越 多 ， 计 算 所 需 的 耗 时 明显 增多 。 在 网 格 划分 为 2x2 的 情况 下 ， 六 种 极 角 分 别 计算 的 


平均 耗 时 为 1.3x103s; 


~ 


网 格 划 分 增加 到 4x4 时 ,六 种 极 角 分 别 计算 的 平均 耗 时 也 随 之 增加 到 2.2x103s; 在 


6x6 的 更 多 网 格 划 分 情况 下 ， 六 种 极 角 分 别 计算 的 平均 耗 时 也 上 升 到 了 2.7x103s; 直至 网 格 划分 增加 到 8x8 
时 ， 六 种 极 角 分 别 计算 的 平均 耗 时 为 7.4x103s。 图 8 中 可 知 ， 网 格 数目 从 6x6 增加 到 8x8 时 ， 耗 时 的 变化 
尤为 突出 ， 六 种 不 同 的 极 角 的 耗 时 都 会 增加 至 少 4.6x103s。 


图 9 展示 了 不 同 极 角 在 不 同 网 格 数目 下 计算 所 得 的 误差 。 不 论 选择 何 
精度 都 有 显著 增长 。 其 中 GAUS 求 积 组 
2 比 3 更 加 精确 


DRZ, ii 
选择 极 角 数目 
5 内 发 生变 化 。 但 由 


Y LÍ 
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S 
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图 9 不 同 网 格 和 极 角 下 的 计算 误差 


Fig.9 Calculation error under different grids and polar angles 


图 6 和 


3.2.2 方位 角 


显 


耗 时 有 
目 


的 增长 ， 


化 。 这 是 由 


方位 角 


Azimuth angles 


8 
12 


16 


为 8 的 3 倍 。 在 内 存 开销 
差 可 知 ， 随 着 方位 角 数 目的 增多 ， 计 和 售 
于 DRAGON 5.0.8 程序 ! 


在 基于 网 格 划 分 6x6， 极 角 选 择 GAUS(2) 
分 析 。 不 同方 位 角 下 的 内 存 开销 


但 到 了 方 


o 在 网 格 划分 分 别 为 6x6 与 8x8 时 , 364 
图 7 可 知 ， 此 时 网 格 划分 8x8 


、 耗 时 以 及 计 香 
ff y 30 之 后 耗 时 几乎 不 再 
中 可 得 出 ,改变 方位 角 的 数目 内 存 开 和 
精度 也 会 愈 发 精确 ， 
方位 角 选 择 的 上 限 为 3 


要 比 其 他 求 积 组 精确 


顾 计 算 精 度 和 效率 。 


— 


种 极 角 类 型 和 个 数 ， 随 着 网 格 划 
至 少 5x104, TE GAUS 求 积 组 中 


FE GAUSC) 极 角 的 计算 精度 仅 在 5x10- 
的 内 存 开销 是 网 格 划 分 6x6 的 1.6 倍 ， 而 耗 时 更 是 
网 格 划 分 6x6 的 2.7 倍 。 根 据 不 同 网 格 划 分 和 极 角 选择 下 对 计算 精度 、 计 算 效 率 和 内 存 开 销 的 分 析 ， 
网 格 划分 6x6， 极 角 选 择 GAUS(2)， 能 较 好 的 


4E 


导出 


的 参数 条 件 下 ， 通 过 不 同 的 方位 角 个 数 完成 进一步 的 敏感 性 
误差 如 表 3 所 示 。 从 耗 时 可 得 出 ， 随 着 方位 角 数 目的 增多 ， 
增长 ， 且 方位 角 数 目 为 30 时 耗 时 约 为 方位 角 数 
省 的 变化 在 5 多 以 内 。 根 据 表 
晶 在 方位 角 数 目 为 30 之 后 ， 精 度 便 不 再 发 生变 


3 中 的 计算 误 


表 3 不 同方 位 角 下 的 计算 结果 


0， 故 综合 考虑 方位 角 的 最 优 参数 为 30。 


Table3 Calculation results under different azimuth angles 


耗 时 
Elapsed time / s 


890.06 


1345.66 


1729.97 


内 存 开销 
Memory overhead / kB 


268492 


271336 


281766 
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相对 误差 


Relative error / 10? 


51.16 


21.83 


13.23 
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30 2671.79 271202 
32 2723.28 271056 
64 2537.84 271412 


3.2.3 特征 线 密 度 


5.90 


5.90 


5.90 


以 网 格 划分 6x6， 极 角 选 择 GAUSO), 方位 角 为 30 的 情况 下 ， 选 择 不 同 的 特征 线 密度 进行 计算 ， 得 到 


因此 随 着 特征 线 密度 的 增 


cn 


hr 


肖 和 计算 时 间 。 


内 存 开 销 、 耗 时 及 计算 误差 如 表 4 所 示 。 不 同 特征 线 密度 下 的 计算 精度 者 
线 密 度 为 12 Cem 的 精度 最 好 , 误差 仅 为 5.9x10-5。 同时 由 于 DRAGON 程序 采用 了 特征 线段 合并 技术 ， 

加 ， 内 存 不 会 出 现 明显 的 增加 R834。 虽然 合并 了 部 分 特征 线段 ， 但 仍 会 使 得 整体 过 
旦 规模 增 大 ， 从 而 导致 耗 时 增加 。 因 此 ， 特 征 线 密度 选择 12 Cem 能 使 计算 误差 最 小 ， 同 时 兼顾 内 存 开 


BHAR, EE 3x104 以内。 其 


表 4 不 同 特征 线 密度 下 的 计算 结果 


Table4 Calculation results are obtained under different characteristic linear densities 


特征 线 密度 


Characteristic linear densi 


8 

10 
12 
15 
20 


3.2.4 功率 分 布 


耗 时 内 存 占用 相对 误差 
ties / cm! Elapsed time/s Memory overhead / kB Relative error / 10? 
1807.41 271382 14.16 
2409.81 272808 10.87 
2662.86 271268 5.90 
2856.16 273492 22.50 
3487.05 282894 12.30 


特征 线 密度 12 Ccm!) 作 


基于 OpenMC 得 到 的 归 一 化 功率 分 布 如 图 10 所 示 。 同时 以 网 格 划 分 6x6、 极 角 GAUSO 方位 角 30、 


为 选 定 参数 ， 计 算 所 得 的 功率 归 一 化 分 布 如 图 11 所 示 。 由 图 11 可 知 ， 计 算 结 果 


与 参考 解 吻 合 较 好 ， 其 ， 
列 ， 与 参考 解 的 相对 误差 
其 值 为 1.61x105, 与 参考 角 


最 大 归 一 化 功率 为 2.506， 位 于 左上 和 角 UO» 燃料 组 件 中 的 第 4 fT 5 列 和 第 5 行 4 
为 0.04%。 位 于 左下 和 角 UO: 燃料 组 件 的 裂变 室 第 26 fT 26 列 的 归 一 化 功率 最 小 ， 


坚 的 的 相对 误差 为 0.32%。 通 过 式 (15 ) 计 算得 到 整个 基 Y 


仁 题 的 均 方 根 误 差 为 0.00255， 


由 此 可 知 ， 该 参数 的 计生 


结果 与 参考 解 吻 合 良 好 。 
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15 
1.0 
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图 10 OpenMC 栅 元 归 一 化 分 布 
Fig.10 Normalized cell distribution in OpenMC 


2.5 
2.0 
15 
1.0 
0.5 
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图 11 DRAGON 栅 元 归 一 化 分 布 
Fig.11 Normalized cell distribution in DRAGON 


通过 改变 不 同 参数 进行 敏感 性 分 析 可 得 ， 在 计算 精度 上 : 影响 最 大 的 因素 是 极 角 的 选择 ， 其 中 GAUS 
求 积 组 要 比 其 他 求 积 组 精确 至 少 Sx104， 由 公式 〈8) 可知， 因为 不 同 求 积 组 会 使 得 角度 离散 过 程 中 的 极 角 
求 积 系 数 w, 发生 变化， 进而 影响 最 终 的 计算 结果 R99。 在 计算 成 本 上 : 网 格 划 分 对 计算 内 存 开销 以 及 计算 时 
间 的 影响 最 大 , 原因 在 于 增多 单个 栅 元 的 网 格 数目 , 会 使 整个 51x51 的 组 件 中 每 个 栅 元 的 网 格 划 分 都 增多 ， 
从 而 显著 增加 其 计算 规模 ， 最 终 影响 到 计算 所 需 的 总 时 间 以 及 内 存 开销 。 


4 结语 


本 文 在 分 析 比 较 MOC、CP 和 IC 计算 原理 的 基础 上 ， 进一步 基 于 2D C5G7-MOX 基准 题 进行 了 数值 计 
算 分 析 , 计算 结果 表明 : 特征 线 方法 在 计算 精度 、 计 算 效率 和 内 存 开销 方面 均 优 于 碰撞 概率 法 和 界面 流 法 。 
其 中 特征 线 方法 的 耗 时 为 23.9min， 只 占 界 面 流 法 的 /16。 特 征 线 方法 内 存 开销 约 37.3SMB， 仅 为 界面 流 法 
的 1/33。 而 碰撞 概率 法 的 耗 时 、 内 存 开销 则 远大 于 以 上 两 种 方法 ， 但 特征 线 方法 和 碰撞 概率 法 的 与 基准 值 
的 相对 误差 都 约 为 6.1x104， 远 小 于 界面 流 法 的 相对 误差 。 因 此 ，MOC 在 计算 2D C5G7-MOX 这 类 强 非 均 
匀 性 组 件 问题 时 ， 其 计算 精度 、 计 算 效 率 和 内 存 开 销 总 体 优 于 碰撞 概率 法 和 界面 流 法 。 

进一步 对 特征 线 方法 进行 初步 参数 敏感 性 分 析 ， 选 择 网 格 划 分 、 方 位 角 、 极 角 类 型 与 个 数 、 特 征 线 密 
度 进 行 计算 分 析 。 通 过 不 同 参数 组 合 下 的 计算 结果 比较 了 计算 精度 、 计 算 效 率 和 内 存 开 销 。 最 终 得 出 网 格 
划分 6x6、 极 角 GAUSQ). 方位 角 30、 特 征 线 密 度 12cm- 是 一 组 能 够 较 好 的 兼顾 计算 精度 和 计算 成 本 的 参 
数 ， 该 参数 组 合 下 的 计算 耗 时 为 45.4min， 内 存 开销 264.7MB ， 相 对 误差 为 5.9x105， 归 一 化 后 的 栅 元 均 方 
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根 误差 为 0.00255， 与 参考 解 吻 合 恨 好 。 由 于 参数 的 选择 与 计算 精度 、 计 算 效 率 和 内 存 开销 间 的 关系 复杂 ， 
各 自 权 重 难以 精确 量化 评估 ， 后 续 研 究 中 拟 引 入 深度 学 习 算法 对 参数 选择 及 优化 进行 深入 研究 。 


作者 贡献 声明 ”到 瀚 源 负责 建立 模型 ， 研 究 并 分 析 数 据 以 及 起 草 论文 ， 张 彬 航 负 责 设计 具体 研究 内 容 和 
方向 ， 提 出 合理 化 研究 方案 以 及 提供 理论 文 撑 ， 同 时 指导 程序 的 使 用 以 及 数据 的 采集 ， 对 论文 的 知识 性 内 
容 做 审阅 并 且 负 责 对 论文 的 修改 ， 袁 显 宝 负责 过 程 监督 以 及 成 果 验 收 ， 张 永 红 负 责 论文 审阅 ， 唐 海 波 负责 
技术 支持 。 
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